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Resizable elastomeric nanopores have been used to measure pulses of ionic current caused by car- 
boxylated polystyrene spheres of diameter 200 nm and 800 nm. The nanopores represent a novel 
technology which enables nanoscale resizing of a pore by macroscopic actuation of an elastomeric 
membrane. Three different pores were employed with variable applied strain, transmembrane poten- 
tial, particle concentration and sphere radius. Theory describing current pulse magnitude has been 
extended to conical pore geometry. A consistent method for interpretation of data close to the noise 
threshold has been introduced, and experimental data has been used to compare several methods for 
efficient, non-destructive calculation of pore dimensions. The most effective models emphasize the 
absolute pulse size, which is predominantly determined by the opening radius at the narrowest part 
of the roughly conical pores, rather than the profile along the entire pore length. Experiments were 
carried out in a regime for which both electro-osmotic and electrophoretic transport are significant. 
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1 Introduction 



Resizable elastomeric nanopores represent one of the most interesting new technologies in the burgeoning 
field of nanopore science. Individual nanopores in thin membranes are attracting interest due to their 
application to sensing of single molecules or small particles, their simplicity, and the growing capability 
in related fabrication and characterisation techniques. Resizable nanopores are each fabricated in an 
elastomeric membrane, so that the nanoscale dimensions of the pore can be altered by stretching and 
relaxing the membrane on macroscopic scales. An initial study demonstrated that individual double- 
stranded DNA molecules could be gated using this mechanism [l|. More recent work has included 
systematic ionic current measurements and analysis of pore actuation 2 , 0] , efforts to characterize the 
pores using AFM and SEM 2, 4], and an initial compilation of theory, experiments and ideas relating 
to resizable nanopores [5]. Resizable nanopore technology is being uniquely developed by Izon Science 
(Christchurch, New Zealand), who supplied specimens, analysis software and the qNano™ actuation 
platform for the present work. 

The fundamental, nanoscale functionality of resizable nanopores lends itself to a wide range of po- 
tential applications. Entirely new processes could arise from this concept, such as mechanical trapping 
of small particles or molecules, controlled mechanical tuning at the nanoscale, or localisation and con- 
finement of reaction chemistry. However, the most immediate interest has been generated by the new 



functionality added to app 



ications which have been studied using conventional (static) pores, such as 



translocation of particles 6l4ldl|. Work towards fast genomic sequencing, perhaps the most exciting po- 



ll 



16 ], has arisen from translocation experiments using nucleic acids 



tential application of nanopores 

With a resizable pore, translocation rales can be con I colled, particles can he gated, and the physics 
of translocation can be studied in novel ways. For translocation measurements, nanopores are filled 
and surrounded by an aqueous electrolyte, allowing measurement of electric current when a potential is 
applied across the membrane. 

Resizable nanopores have further practical advantages over static pores, which have employed a 



range of membrane materials and nanopore fabrication methods 



18 



191 ]. The most widely-used 

conventional pores are 'solid-state' nanopores fabricated in thin, rigid, silicon-based jf| or polymer [3, 0] 



2 



14 



15 



membranes and 'biological' pores {20j , especially the a-haemolysin pore 
membrane [g, 



21] supported by a phospholipid 



22l | . In comparison with such pores, resizable nanopores are less likely to irreversibly 
clog, as the elastomer can be stretched to allow the pore to clear. The elastomer is light and the effective 
pore size for a single specimen can be varied over approximately one order of magnitude |2], enabling 
it to perform the tasks of several static pores. The fabrication method is quick, simple and relatively 



inexpensive 



Elastomeric pores are chemically, mechanically and thermally stable in comparison 
with biological pores, which must be carefully supported under laboratory conditions. The major present 
drawback of solid state pores, including resizable nanopores, is poor reproducibility of geometry below 
~50 nm diameter 



231 ] . It is hoped that studies such as the present one will pave the way for the 
manufacture, characterization and widespread use of reliable, molecular-scale resizable pores. 

In this article, we present experiments in which polystyrene nanospheres have been detected using 
resizable nanopores. The work demonstrates that this technology can be used for translocation experi- 
ments, in a manner similar to static nanopores. The experiments build on initial work [5] by obtaining 
data using various nanopore specimens, applied strains, applied potentials, particle sizes and particle 
concentrations. The investigation is the first to consider the relationship between current signals, noise, 
particle size, particle type and pore geometry with relation to resizable pores. We address the fact that 
particle-induced current signals may not always indicate passage of the particle from one side of the 
membrane to the other. An important feature of this work is a direct comparison of existing theoreti- 
cal models for relating the resistance pulse size to pore dimensions P, I23U28]. including a conical-pore 
approach developed in the theoretical background. An efficient, non-destructive model for determining 
pore size would be valuable for many envisaged applications of this technology. Particle transport is also 
specifically addressed, both in experiments and the identification of relevant aspects of the Space-Charge 



model (e.g. 



24, 



291 ]) in the theoretical background. 
Resizable nanopore specimens are currently available for sensing particles in the hundreds of nanome- 
tres range; the size of the polystyrene particles used here (200 and 800 nm diameter) indicates the size 
of the narrowest part of the pore. This size range is relevant to many viruses of interest in human and 
veterinary medicine, agriculture and environmental studies. Many of the recent developments in virus 
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detection technology have required fluorescent molecular tagging of the virus. The most relevant 



ous work includes recent nanopore studies using particles of similar size and material 



23, 



26, 



31 



33, 



25, 



27 



previ- 



3JH3J] as 



34], which included 



well as earlier work relating to the development of Coulter counters 
detection of particles as small as 60 nm [34]. The membranes used in experiments are relatively thick 
(200 fim when unstretched) in comparison with other studies. Pores are conical 2, 4, j|, and some 



current rectification, as studied in some depth elsewhere 



35H40H , has been observed previously [2] . The 



pores used do not close entirely when the membrane is relaxed. 

2 Theoretical Background 

In this Section, theory relating to two major topics is presented. Firstly, we recount descriptions of 
transport and apply them in the context of the present work. Th e approaches used have been widely 



applied to both ions 



30, 



31 



41 



13] and pari iclcs ^, 



10 



23 



25 



30 



311 ] in e: 



elsewhere. Secondly, we present theory developed by deBlois and Bean 



studies y,|23j-l25, 



27 



xperimental nanopore work 



261 ] and widely used in recent 



281 ] to describe the size of a resistive pulse generated by a spherical particle passing 



through a nanopore. This theory is reviewed, so that assumptions can be tested by the experimental 
work, and we explicitly extend the existing theory to conical geometry. 



2.1 Transport of Charged Particles 



Several papers 



24, 



29, 



30, 



38, 



43] summarise theoretical approaches that have been used to study 



transport of ions and particles in pores, ranging from biological ion channels to micron-sized synthetic 



pores. We employ a simple approach 4J] devised for assessing the relative contribution of various 
transport mechanisms, and their functional dependencies, for experimental work. This approach is a 



24 



29]. 



simplification of the comprehensive, but computationally expensive Space-Charge model 

Transport of aqueous ions or larger charged particles is described by the Nernst-Planck equation. The 
particle flux J is given in terms of diffusive Jdiff > electrophoretic J ep and convective J C onv contributions, 

by 



4 



J — Jdiff 4~ Jep + Jconv 

= -DVC + -——DC~E + Cv. 

Here, I? is the diffusion coefficient, C is concentration, ks is Boltzmann's constant, T is temperature, E is 
the electric field and v is the convective flow velocity. The electronic charge magnitude e is multiplied by 
£ to find the total effective charge on the ion or particle; £ = — 1 for an unscreened electron. Alternative 
notations use the gas constant R = IzbNa and the Faraday constant F = N^e, where Na is Avogadro's 
number. At 293 K, a dimensionless ratio for determining the significance of the diffusion term is [441 ] . 



'diff 



O025 

Jep ^ { ' 

Generally, it is expected that |£| > 100 for carboxylated polystyrene nanospheres of 0(100) nm diameter, 
in which case the ratio in Eq. [2] is < 5% for any potential above ~ 0.01 V. Therefore, diffusion should 
not significantly contribute to particle transport in the present work. Electrolyte concentration is the 
same on both sides of the pore, so ionic diffusion is also not important. 

In the absence of other external forces, such as a pressure gradient, electro-osmosis is the sole 
mechanism for convective transport. Referring to the geometry defined in Fig. [H and using the simplified 
model of a l ong , thin cylindrical pore (o = b = do, I » ao) to model the electro-osmotic flux (J z = |J|), 



we find 



44 



45J 



Jz.e.o — ; L'-G/z 1 



Attt] \ kclI («o) y , s 

Here, r\ is the dynamic viscosity of the fluid, e is the dielectric constant, Vo is the potential of the pore 
wall and k _1 is the characteristic thickness of the electrical double layer at the pore wall, and I n is the 
nth-order modified Bessel function of the first kind. The subscript eo refers to electro-osmotic flow. The 
value of A increases asymptotically towards 1, and is greater than 0.9 for a > 20 nm when using 0.1 M 
KC1 = 9.62 x 10 -8 m _1 [46]). The ratio of electro-osmotic to electrophoretic particle fluxes is 



b f 




FIG. 1. A schematic cross-section through a truncated circular conical pore, showing geometrical quan- 
tities used in the theoretical analyses. Cone entry radii are a and b and the pore length (equal to the 
membrane thickness) is /. A translocating sphere of radius a' is shown in the conical passage and a 
cylindrical polar co-ordinate system is defined. 



J, 



3a'etpo 



A, 



(4) 



'z,ep 2£e 

where the particle is assumed to be a sphere of hydrodynamic radius a'. For water at 293 K, this ratio 
is —0.10 A, when calculated using values of a' = 100 nm, = — 75 mV and £ = —500. This calculation 
demonstrates that both transport mechanisms are important in the present experiments, while suggesting 
strong dependence on surface charges and potentials, which are difficult to experimentally determine. 
The value tpo = —75 mV is equivalent to surface charge density of 0.3e nm~ 2 in 0.1 M KC1, consistent 



with literature values for track-etched polymers 



35 



37 



33). Note that 



Jz,ep 



CE Z 



(5) 



67rr]a' Atttj 

so both contributions are linearly dependent on particle concentration and potential difference. 

A corresponding comparison of electro-osmotic and electrophoretic contributions to ionic current 
carried by ions 



44, 



45] reveals that electrophoresis dominates in the experimental regime considered 
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here. If the second (electrophoretic) term dominates in Eq. [TJ pore size can be estimated from ionic 
current measurements by treating the pore as an isotropic, homogeneous conical conductor with bulk 
conductance p equal to that of the electrolyte (see Eqs.E anc [1] below). This approach has been widely 



used 9, 



25, 



26 



28, 



30, 



32], including with conical pores 2j, |7|, |8|, 



Current rectification has been observed using resizable nanopores, and should be studied further, 



althoug 
shown 



i the polarity of applied potential was not varied during the present experiments. Studies have 



3514401] that rectification is determined by asymmetric geometry, non-uniform surface charge 



distributions and ionic concentration gradients. In previous work, the current anisotropy observed using 



elastomeric nanopores 



was at most 25% for a transmembrane potential of +200 mV. 



2.2 Translocation Event Size 

When a particle travels through a nanopore at constant applied potential, a momentary change in 
pore resistance is observed - a resistance 'pulse' caused by a translocation event. The characteristics 
of this pulse yield information on the size, shape and nature of the particle. Qualitatively, it is easy 
to understand that an insulating polystyrene sphere will 'block' the pore, increasing its resistance, as 
it passes through the membrane; that a smaller sphere will produce a smaller increase in resistance; 
and that the resistance increase is greatest when the particle passes through the narrowest part of the 



pore. Nearly forty years ago, deBlois and Bean 



used classical electrostatics to derive quantitative 



expressions for pulse magnitudes that have been widely applied in recent studies. In this section, 
approaches detailed by these researchers and others are recounted, and in some cases extended to 
incorporate a linear conical geometry which is relevant in the present work. It is of interest to compare 
the various models using experimental data in order to probe the functional relation between particle 
size, pore size and pore geometry. 

DeBlois and Bean's approaches are approximations, each of which assumes that the particle is spher- 
ical and that the surfaces of the membrane, pore and translocating sphere are uncharged. The full 
solution of Laplace's equation, which is required in order to solve the classical electrostatic problem, is 
not even available for the simple geometry of an insulating sphere within a conducting tube (Fig. [1]). 
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Moreover, a real system introduces additional complications such as non-cylindrical pore geometry, 
solution and surface chemistry effects and aspherical particles. 



2.2.1 Simple Electrostatic Approaches 



Recent studies of spherical proteins [9|, 



28] and other particles 9|, 23II25I. have based their interpreta- 



tion of translocation pulses on the simplest of the approaches presented in [26[ , which assumes that the 
sphere radius a' is small in comparison with the cylindrical pore radius a. For the geometry defined in 
Fig. [H we initially consider a pore which is cylindrical (a = b = ao) and very long (I » a). In this case, 
the electrical resistance of a pore filled with electrolyte of resistivity p is 



R, 



cyll 



pi 



■nan 



(6) 



Maxwell's expression for the resistivity of a solution containing insulating spheres p e ff with volume 
fraction / is 



Peff = P 1 



3/ 



(7) 



where / is just the ratio of the single particle volume to the pore volume, 



fcyl 



4a /3 
3op' 



(8) 



The absolute change in resistance during a translocation event, or the amplitude of the resistance pulse, 
is 



&R C yl — Rcyl2 ~ Rcyll 



2 pa 



/3 



(9) 



where R cy i2 is the equivalent resistance when the pore contains one spherical particle. Dividing Eq. [9] 
by Eq. [6l the fractional change in resistance is 



AR cyi 2a 



13 



Rcyll CLqI 



(10) 
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For spheres of constant size, the absolute resistance change varies as and the fractional resistance 
change varies as 2 . Therefore, current pulses are most clearly observed when the pore size matches 
particle size as closely as possible. 

The same approach can be extended to different pore and particle geometries, while still assuming 
that ao » a'. If I » ao does not hold, end effects are significant, in which case the solution given in 
Eqs. E]to[TU]is simply modified by replacing the pore length I with the approximate factor I + 1.6ao [26]. 
For the present work, we explicitly consider a linear conical pore with end radii a and b, in which case 
the pore resistance is 



the volume fraction is 



R, 



conl 



ixab^ 



4a' 3 



(11) 



(a 2 + ab + b 2 )V 



(12) 



and 



6pa> 

£XKcon - nabi^ + ab + b 2 )' [ 6) 

Equations [11] to [13] are derived using volume integrals over a conical pore, and apply when a, b and a' 
are much smaller than I. For a truncated cone in which end effects are significant, note that the pore 
resistance is simply the series sum of the contributions from the ends and central part of the pore. By 
analogy with the approximation for the cylindrical pore, we can replace I with / + 0.8 (a + b) in Eqs. [11] 
and [12] We obtain a conical version of Eq. [10] by dividing Eq. [13] by Eq. [TTl 



AR r 



6a' 3 



Rconi I {a 2 +ab + b 2 ) ' 



(14) 



2.2.2 Other Approaches 



Ot 



m 



rer methods may be more applicable in the large-sphere cylindrical-pore limit. An approach presented 
26J and attributed to Gregg and Steidley 



481 ] has been used in some experiments with latex spheres 



{27]]. In this method, pore resistance is calculated by integrating the annular section surrounding a 
spherical particle over the length of the pore. This approach does not require that the particle is 
significantly smaller than the pore size, but necessarily underestimates the pore or blockage resistance, 
because a uniform current distribution is assumed, whereas any nonuniformity creates a larger resistance. 
The absolute change in pore resistance is 



Aft, 



cyl 



2p_ 
vra 



sin 



(ao) 



1 

2\ 2 



a 
a 



(15) 



In the limit ao » a' , this result differs from the previous small sphere result (Eq. E|) by a factor of 



26]. 



two-thirds. However, Eq. [15] is asymptotically valid when the sphere size approaches the pore size 

DeBlois and Bean more rigorously approximated the problem using a different method, in which 
the calculated electric field bulges around the insulating sphere contained within the pore. Deflection 
of field lines is minimal, especially for a small sphere, allowing calculation of an upper limit for AR cy i. 
This method applies for any sphere size when I » ao, and gives a more rigorous result than Eq. [TBI in 
the regime a' « ao, asymptotically approaching the result from Eq. in the small sphere limit. The 
result may be most appropriate for the intermediate region in which the sphere is smaller than, but of 
comparable dimensions to, the pore size. It can be approximated using the series expansion 26] 



AR 



cyl 



2 pa 



/3 



Zimmermann and Jelsch 



,/3 



,/6 



1 + 1.26- 



1.1- 



(16) 



33j ] considered resistance pulses generated by cells, and therefore modelled 
aspherical particles of finite conductivity. Using their approach, the fractional volume change for a single 
particle in a cylinder is 



AR 



cyl 



R, 



■cyll 



fsft 



4a 



/3 



3oq/ ' 



(17) 



where f s and f c are shape and charge factors respectively. This result reduces to Eq. [UJ for a perfectly 



insulating spherical particle, when f s = | and f c = 1. The value f s = 1 applies to an extended 
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cylindrical particle, similar to an elongated polymer. For a perfectly conducting sphere, the charge 
factor is f c = — (f s — 1) so there is a continuum of possible fractional resistance amplitudes which 
includes zero. This analysis gives some simple insight into the dependence of translocation current peaks 
on charge distributions, and the associated solution chemistry. Polystyrene is an effective insulator, but 
surface charges, polarization charges and associated screening charges are present in any nanofluidic 
system. For example, screening can reduce the value of £ by more than two orders of magnitude for 
carboxylated polystyrene spheres [49]. 



3 Experimental 

Experiments were carried out using a 0.1 M KC1 solution, prepared using deionised water (18.2 MO) 
and buffered at pH 8.0 using 0.01 M tris base (pK a 8.1). Small amounts of Triton X-100 and EGTA 
were added to aid wetting and as a chelating agent respectively. Reagents were purchased from BDH 
and Sigma. Particle concentrations are calculated relative to 'undiluted' (as-received) aqueous solutions 
of carboxylated polystyrene spheres (1% wt/v, Corpuscular). Nominal particle diameters of 200 nm and 
800 nm were used. 

Three resizable pore specimens were used. Each pore is located at the centre of a thermoplastic 

i' q, y, b (Fig. 



polyurethane (TPU) 'cruciform' 



2(a)). They are fabricated by puncturing the 200 ^m- 



thick central septum using a specially etched, sharpened tungsten tip. The geometry of each pore is 
roughly a truncated circular cone, tapering from a large hole on the 'cis' surface to a smaller aperture 
near the 'trans' surface. For pores manufactured under the same conditions as the specimens used here, 
SEM images have revealed cis surface openings of characteristic radius 15 //m 2j, [4], [5J. AFM imaging 
[2, 5] suggests that the narrowest part of the pore, where the electronic sensing is most sensitive, is 
within ~1 fim of the trans surface opening, rather than on the surface itself. 



Pores are actuated by extending the dista. 
to stretch radially, with azimuthal symmetry 



ends of the cruciform legs, causing the central membrane 
2J. Each specimen undergoes stress-softening prior to use, 
so actuation is largely reversible and reproducible outside the near-pore region, which is catastrophically 
overstretched during pore formation 5J. Following actuation, the cis side pore radius, and the thickness 
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FIG. 2. Apparatus used in the experiments. The circular central septum of a TPU cruciform (Fig. 2(a) 



with a millimetre scale), where nanopores are located, is 200 /im thick. The remainder of the cruciform 
is predominantly ~0.8 mm thick. Cruciforms are placed within a fluid cell (Fig. |2(b)[ ) so that o-rings 
form a loose seal on the cruciform surface. The lower half of the cell is enclosed, whereas solutions can 
be added or 'flushed' into the open upper surface. The cell is placed on the top of the Izon qNano 



(Fig. 2(c) ) so that teeth fit in the holes at the ends of the cruciform legs. These teeth are actuated by 
turning the dial on the side of the qNano. 
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of the membrane, can be estimated by assuming that the material is Hookean, elastic and incompressible 
Q]. Improved assessment of the relationship between macroscopic strain and deformation local to the 
pore, including deviations from a linear stress-strain profile, is the subject of ongoing work. 

A qNano system (Izon Science, Christchurch, Fig. [2]) was used to actuate the cruciforms and carry 
out measurements. The central part of each cruciform was wet by electrolyte using a fluid cell which 
allowed loading of fluid into a reservoir below the specimen, while solutions containing particles could be 
pipetted on to the open top half of the cell. Vernier callipers were used to measure the distance between 
the outer rims of holes on opposite cruciform legs. The cruciform 'stretch' (see Table IJ) is defined as 
the difference between this measurement in the stretched and relaxed states, where the latter value is 
42 mm. Current was measured and voltage applied using Ag/AgCl electrodes in contact with each half 
of the fluid cell. The current signal (e.g. Figs. 3(a) and |3(b) | was digitized at 10 kHz using the qNano's 
built-in electronics, while software customised for the qNano (Izon Instrument Control System v2.1) was 
used to set the transmembrane potential, identify events from the current signal and collate the recorded 
data into histograms. 

Typical records of raw data and examples of event size distributions are shown in Fig. [3) The 
current and translocation traces are comparable to other work, using streptavadin for example j^j. 
The noise threshold for event size, below which events are ignored, is set at 0.05 nA, consistent with 
the characteristic peak-to-peak baseline signal noise. This current threshold is higher than in other 
nanopore studies in which molecular-scale particles have been detected. We emphasize that these are 
initial experiments using comparatively large pores with a portable, desktop system under ambient 
laboratory conditions. In control experiments with no particles present in the electrolyte, a maximum 
of 1 event s _1 was identified by the qNano system. All such events had magnitude between 0.05 and 
0.06 nA. When nanospheres are introduced to one half of the fluid cell, increased noise is generated 
by particle activity, such as partial blockages of the pore and interactions of spheres with each other 
and the pore walls. Fig. [3] demonstrates that events near the noise threshold may be more frequent 
than a histogram peak caused by translocations at greater event magnitude, or less frequent, as in 
Figs. 3(c) and |3(d)] respectively This translocation peak can be dominated by noise (as in Fig. 3(c)), 
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Time / s 



Time / s 



(a) 



0.4 0.6 

Event Size / nA 



0.4 0.6 
Event Size / nA 



(c) 



(d) 



(b) 





0.4 0.6 
Event Size / nA 



(e) 



FIG. 3. Examples of raw data and event histograms. Figures 3(a) and |3(bj| show excerpts from typical 
raw current records using cruciform B (see Table [I]). Typical events, indicated by a dashed ring on 
the main trace, are expanded in the inset traces. In Fig. |3(a) from series 1, events were caused by 



800 nm spheres and the applied potential was 0.42 V; Fig. 3(b) is from series 2, using 200 nm spheres 



and the potential was 0.44 V. Figures |3(cj| [3(d)1 and [3(e)] are event size distribution histograms from, 
respectively, cruciform C series 1 at 0.40 V (1 minute duration), cruciform B series 1 at 0.50 V and 
cruciform B series 1 at 0.30 V (both 5 minutes duration). Event sizes were characterised relative to the 
smoothed background current in the vicinity of each event [5J and collated into 10 pA bins. 
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especially for smaller particles, due to the dependence of event size on particle volume in the simple 
model (Section I2.2p . In such cases, the tally of events is taken to represent particle 'activity', rather 



50] 



than just translocations; this approach is consistent with similar work elsewhere 

The two key parameters extracted from the raw data for the Results are (i) the size of translocation 
blockage signals and (ii) the frequency of events. The histograms demonstrate that it is not useful to 
calculate the translocation event magnitude using a mean of all data. It is difficult to separate translo- 
cation events from noise events when the tails of the distributions overlap. Therefore, the characteristic 
size of a translocation current blockage is defined here as the modal peak once the peak at the noise 



level is discarded. When there is no discernable translocation peak (e.g. Fig. 3(e) ), the data is usually 
discarded, but the peak at the noise threshold can be interpreted as an upper bound on the actual 
translocation blockage magnitude. This approach is comparable to the identification of 'clusters' used 
elsewhere 9|. 

In order to ensure consistency in frequency data, all identified events above the 0.05 nA noise limit 
are included. Therefore, frequency data represents all activity, not just translocations. Use of this 
approach is supported by the smooth linear relationship between event frequency and concentration 
observed in the Results. The smoothed baseline current adjacent to each event was also recorded in 
order to calculate the fractional event size, and for estimation of pore size using the bulk conductivity 
method. 

Event duration can be determined by considering the current signal around the event in detail, as in 
the inset to Figs. 3(a) and |3(b)j This parameter is not analysed in the present work due to inconsistency 
of the data. As discussed above, it is difficult to discriminate translocations from other events. The 
beginning and end of an event can also be ill-defined, especially when it is considered that the spheres are 
passing through a long cone. Further work is required to address these issues, or to set up experiments 
so that pulse durations are more regular. 

Experiments are grouped in specific 'series' for each cruciform, with each series corresponding to a 
stable cruciform setting at which translocations were observed, as shown in Table [H Within each series, 
one variable, usually voltage, was altered while a number of separate data files ('runs') were recorded 
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Table I. Summary of cruciform specimens, experimental 'series' and the associated experimental condi- 
tions. As described in the main text, concentration is relative to as-received solutions and the cruciform 
stretch is measured relative to a resting length of 42 mm, which is also used to calculate strain. 



Cruciform 


Series 


Stretch 


Strain 


Concentration 


Sphere Radius 


Voltage Range 






mm 




x 10- 4 


nm 


mV 


A 


1 


10.3 


0.25 


2 


400 


0.16-0.30 




2 


10.4 


0.25 


2 


400 


0.16-0.30 




3 


16.0 


0.38 


2 


400 


0.16-0.30 




4 


19.1 


0.45 


2 


400 


0.16-0.30 


B 


1 


19.2 


0.46 


1 


400 


0.34-0.50 




2 


19.2 


0.46 


1 


100 


0.38-0.46 




3 


10.0 


0.24 


2 


400 


0.12-0.50 




4 


10.0 


0.24 


2 


100 


0.16-0.49 




5 


9.51 


0.23 


Varied 


400 


0.30 only 


C 


1 


13.5 


0.32 


2 


400 


0.12-0.46 




2 


13.5 


0.32 


2 


400 


0.14-0.44 



for analysis. For cruciforms A and B (series 1-4), each data point at a particular voltage represents a 
run of typical duration between 5 and 10 minutes. For cruciform B series 5 and cruciform C, each data 
point represents the average value from several runs (typically four or five) , each of typical duration 30 s 
to 1 minute. Voltage was generally varied in random chronological order to avoid systematic error. 
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FIG. 4. Figure 4(a) is a plot of event frequency as a function of concentration (relative to as-received 
colloidal solutions) at 0.30 V applied transmembrane potential (cruciform B, series 5). The two experi- 
ments represent data captured on consecutive days. R 2 = 0.973 and 0.997 for linear fits to experiments 
1 and 2 respectively. Figure |4(b)| shows event frequency plotted over a working range of voltages for 
two series using cruciform C. R 2 values for a linear fit were less than 0.5 for these series, which were 
recorded at constant particle concentration and pore resistance. Plotted points and error bars represent 
the mean and standard error in the mean for multiple recordings. 
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4 Results 



4.1 Dependence of Event Frequency on Concentration and Voltage 



The data in Fig. 4(a) are consistent with the predicted electro-osmotic and electrophoretic particle flux 
(Eq. [5]), as well as previous experiments [5|. It is notable that the trend does not pass through the 
origin, suggesting that the simple transport model of Section \2. II does not extend to low concentrations. 
Data plotted in Fig. |4(b)[ which are typical of other experiments, suggest a weak positive trend between 
event frequency and applied voltage, especially if the outlying data with relatively large error bars is 
discounted. The transport model suggested in Section 12.11 predicts that the event frequency should 
increase linearly with applied voltage. There is little evidence for such a trend; any linear positive trend 
would appear to intercept above the origin. In particular, the significant scatter in the data appears 



inconsistent with the result from Fig. 4(a) 



There are several possible reasons for the inconsistency between Figs. 4(a) and |4(b)| In a broad 



sense, it is most likely that the data reflects differences between the two specimens used. One other likely 
explanation is that pore morphology changes unpredictably with applied voltage, producing the marked 
difference in the scatter of the data in the two figures. The polymeric internal pore surface and adjoining 
electrolyte include regions of aggregated charge that could experience a force in an applied electric field. 
Parts of the pore surface near the constriction can be ragged and mechanically unconstrained due to the 
fabrication method. 

Other factors could explain the nonlinear relation or generally poor data consistency observed in 
Fig. |4(b)[ but not necessarily both. The simplifying assumptions required in order to arrive at Eq. [5] 
should not explain the scatter in the data. Relevant issues would include any electric field leakage near 
the tip of the cone, the suitability of the Space-Charge approach for a conical geometry, and possible 
variations in particle concentration when the pore diameter approaches the particle size. Three further 
experimental factors could be significant. Firstly, the number of recorded events represents all activity, 
rather than just translocations (discussed above). Secondly, apparently spontaneous changes between 
distinct, stable 'modes' are observed when using this device, evidenced by step changes in baseline 
current and possibly caused by adhesion of particles to the polymer surface. However, switches between 
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FIG. 5. Figure 5(a) plots the modal resistance change of events against baseline nanopore resistance. 
Each data point represents the average of data collected at a particular applied potential. 200 nm 
particles are used in cruciform B series 4, for which the data represent an upper bound on translocation 



resistance pulse magnitude (see Section [3]); 800 nm particles are used for the other data. In Fig. 5(b) 
baseline current is plotted against applied potential for experiments using cruciform C. R 2 = 0.97 for 
linear fits to the data in both cases. 



such modes are typically identified if they are significant. Thirdly, it is possible that a small pressure 
difference is introduced across the membrane as it is sealed to the lower half of the fluid cell. 



4.2 Translocation Resistance Pulse Sizes 
4.2.1 Voltage Dependence 



Figure 5(a) shows the modal size of the resistance pulses associated with translocation events. Overall, 
there is a rough positive trend between the size of a resistance pulse and the baseline resistance. The 
exact relation is inconsistent between different series due to geometrical variations relating different 
pores, stretch states and particle sizes. These observations do not depend strongly on the applied 
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transmembrane potential, consistent with the analyses in Section [221 It is apparent in Fig. 5(a) that the 
range of R\ values is limited for any given series. Further, Fig. [5(b)] demonstrates that baseline resistance 
is typically characterised by an Ohmic I-V plot over the range of interest. Although rectification was not 
specifically studied here, the Ohmic characteristic suggests that any current asymmetry is not significant 
in the interpretation of results. 

Experiments using cruciform B demonstrate smaller resistance pulses for smaller particles in two 
series at the same cruciform strain. This trend is intuitive and qualitatively in line with the theoretical 
approaches. Quantitatively, there is a factor of ~10 difference between the modal resistance peak for 
800 nm particles and the upper bound on that peak for 200 nm spheres, using a cruciform held at the 



same strain (Fig. 5(a) ). This is lower than the difference predicted (e.g. Eq. [9]) if the blockade magnitude 
scales with the sphere volume. Taking into account that the peak for 200 nm spheres is an upper bound, 
the data therefore suggest that improved resolution is required to accurately study particles of such 
disparate sizes using a static pore. A key method of improving the signal-to-noise ratio is to match 
the particle size more closely to the pore size, a process which is directly enabled by resizable nanopore 
technology and which will be the subject of further work. 

4.2.2 Estimates of Pore Radius 

The estimated pore radius derived from experimental data, denoted a es t, is heavily dependent on the 
model used for the calculation. At present, this value should be treated as an indication of effective pore 
size for comparative purposes, rather than necessarily an accurate determination of pore radius. Typical 
data showing the calculated estimated radius function of applied potential are shown in Fig. [6j 

These data are calculated using the modal fractional event size (Eq. [TO]) at each potential. In each 
case, a es t is calculated with the nanopore in a stressed state, with changes of the membrane thickness I 
estimated using the simple approach described in [3|]. This plot further explores the assertion of voltage 
independence and an effectively Ohmic relationship observed in relation to Fig. and consistent with 
the analyses in Section [2.21 

For most of the data series in Table [H results are similar to those plotted in Fig. [6] for cruciform B 
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FIG. 6. Estimated pore radius a es t, calculated from the fractional event size using Eq. [TUl is plotted 
as a function of applied potential. Stretch-induced membrane thinning has been accounted for in the 
calculation, and each data point represents the average of all events recorded at a particular voltage. For 
the cruciform C series, error bars represent the standard error in the mean of several, relatively short 
runs at the same potential. 
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Table II. Comparison of the theoretical approaches to calculating a es t, as applied to two data series 
obtained under similar experimental conditions (see Table [I]) . Thinning of the stretched membrane is 
treated as described in jjj. All of the models assume that a es t « I, and uncertainties are dominated by 
choice of the model used for calculations, as discussed in the text. Estimates refer to the smaller pore 
radius in the conical case, and to the narrowest part of the pore when using the event resistance change. 



Model for a est 


Assumes a' « a es {! 


Cruciform C Series 1 
fim 


Cruciform C Series 2 
/im 


From Baseline Resistance: 








Cylinder (Eq. ED 


Yes 


3.5 


3.1 


Cone (Eq. HTJ 


Yes 


0.6 


0.5 


From Event Resistance Change: 








Cylinder (Eq. © 


Yes 


1.6 


1.2 


Cylinder, fractional (Eq. [TO]) 


Yes 


0.7 


0.4 


Cylinder, lower bound (Eq. [T5|) 


No 


1.5 


1.1 


Cylinder, upper bound (Eq. [JBD 


No 


1.6 


1.2 
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series 3 - consistent over the working range of voltages. For cruciform C series 1, there is an indication 
of a positive trend in estimated pore size with voltage. Such a trend is atypical of the data for all 
series, although the degree of scatter is consistent with the other data sets. Comparing the two data 
series for cruciform C, it is clear that apparent pore dimensions can change between series in which 
experimental conditions are repeated. The error bars for the cruciform C series indicate that the results 
for a particular series, applied potential and sphere size are comparable to the scatter over a range of 
potentials and possible differences between series, as indicated by variations in estimated pore size. 

Various methods of calculating estimated pore size a es t from current measurements are compared 
in Table [II] Three models, using absolute resistance data with Eqs. US and HS1 give very similar 
results. It follows that the assumptions differentiating these approaches are relatively insignificant in 
the pore-sizing process. These assumptions are (i) a' « ao, (h) the electric field is uniformly distributed 
for Eq. [T5] and (iii) the electric field bulges as required for Eq. [THJ In contrast, the method of dealing 
with pore geometry (including end effects) significantly varies between these three models and each of 
the other models. 

DeBlois and Bean's simplest analysis can be applied to both absolute (Eq. ED and fractional (Eq. [TU|) 
resistance measurements. The absolute case, which assumes cylindrical pore geometry and does not 
depend on I or b, gives a significantly smaller pore size than the cylindrical result from baseline resistance, 
because the absolute pulse is dominated by the signal when the sphere is within the narrowest part of 
the conical pore - at the small pore end. Values calculated using Eqs. [15] and [16] use a similar geometric 
approach. Results calculated using the fractional resistive pulse data are two to three times smaller than 
the absolute resistance data. This is expected, because a es t cc (R / AR) 2 for the fractional data and 
a es t cc (1 / AR) 4 for the absolute data; the measured value of R, averaged over the length of the pore, 
is relatively low, whereas the value of AR is mostly dependent on the narrowest part of the pore. 

As expected, baseline resistance data gives the largest value for the cylindrical case, because that 
model assumes consistent resistive losses along the pore length. For the conical model, resistive losses 
are concentrated near the smaller end of the truncated cone. The conical calculation uses estimated 
values of the larger pore radius b, determined using SEM images of typical pores at rest and adjusted 
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a. 



to reflect membrane stretching using the approach described in 

Typically, the standard error in measurement of resistance or modal event size is < 5%. Uncertainties 
in the calculated values from Table [TT] are dominated by geometric factors, most of which have been 
discussed above in comparing the different methods of modelling. To illustrate this point, we note that 
the two further conical models (Eq. [13] and [T4"|) give outlying results (a es t < 10 nm) because of their 
dependence on b. The probable explanation for this observation is that the profile of the pore, rather 
than being a linear truncated cone, is convex (trumpet-shaped) 5j|. It follows that the volumetric ratio in 
Eq. [8]is an underestimate. If this explanation is correct, the calculation of a es t using the conical, baseline 
resistance method should also be an underestimate; this is entirely reasonable given that the smaller 
pore radius should accommodate the translocating particles. It should further be noted that the model 
for membrane thinning and pore resizing uses linear elastic material parameters. The approaches 
which use the experimental value AR have the disadvantage of introducing greater fractional random 
error associated with this differential measurement. However, the greatest measurement uncertainty 
in the present experiments is associated with the larger pore radius b, which is typically ~ 15 + 5 //m. 
Improvement to techniques for imaging these pores is ongoing 4]. Other issues, which are less consistent 
with observed differences between model calculations, could arise from the method of electronic sampling 
and from significant charge effects. 

There would be great practical benefits associated with determining pore size and shape using simple 
current measurements. A reliable method would enable expedient use of the technology, without the need 
for expensive, time-consuming and potentially destructive methods of specimen characterisation. Of the 
approaches considered, the most reliable are likely to be those models using absolute pulse data with 
cylindrical geometry. These models estimate the pore size where the particle is most constricted within 
the pore, so the value of a es t should refer to the characteristic pore size at or near this constriction. 
These approaches do not require detailed knowledge of the pore profile, or the values of b and I, so 
modelling of stretch effects is not required. Improved control and understanding of pore geometry will 
be key to developing pore-sizing techniques, while charge effects are likely to become more significant 
for smaller pores. 



24 



Further development of the analytic current blockade analyses could extend the models for which 
a 1 « ao does not hold to a conical pore, in which case the effective solution resistance is not constant 
down the length of the pore. A range of more realistic, non-conical profiles such as those used by 
Ramirez et al. [38] could be considered. End effects might be important when a particle is near the 
most active, narrow end of the nanopore. The usual assumption for neglecting the pore ends is a « I, so 
significant effects are likely when the particle is within a few microns of the narrow pore opening. The 
total population of spheres within the cone of the pore could also be considered. The presence of multiple 
spheres could significantly impact the effective solution conductivity and current noise characteristics. 
Finally, surface charge effects should be included for application to future experiments using smaller 
resizable pores. 



5 Discussion and Conclusions 

This paper has presented initial experimental work using novel resizable nanopore technology. Vari- 
ation of a number of experimental variables has been demonstrated using theory and experiments. 
Experimental trends agree with those established using other apparatus. For example, the absolute 
size of experimental translocation resistance pulses increases with the baseline resistance, meaning that 
translocation signals are clearest for a particular particle set when the resizable pore radius is as small 
as possible. The increase in pulse size between signals for 200 nm and 800 nm particles is consistent 
with volumetric scaling of pulse size. There is little dependence of pulse size on applied voltage, a result 
consistent with a predominantly Ohmic I-V relationship, which was independently verified. The exper- 
imental data also demonstrated some inconsistency, manifested as scatter, reflecting that measurements 
were carried out using a dynamically variable system. 

A method for accurately estimating nanopore size from current measurements alone would be valu- 
able for applications of this technology. For that purpose, several theoretical models were directly 
compared using experimental data, including a novel approach in which truncated cone geometry was 
applied to resistance pulse magnitudes. The most accurate models rely on absolute resistance pulse 
height, emphasizing data recorded when the sphere is at the narrowest part of the pore. Such models do 
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not heavily depend on the pore profile through the entire membrane thickness. Overall, the spread of 
values calculated using different models demonstrated that internal pore geometry is the most important, 
least well-understood experimental parameter. 

Several issues relating to experimental configuration and data collection have been directly addressed 
in this work. It was established that both electrophoretic and electro-osmotic flow were probably sig- 
nificant for particle transport in the experiments described. Key functional relationships presented in 
the theory section allowed effective handling of this experimental regime. In particular, it was predicted 
that particle flux should be proportional to concentration and electric field. Experiments supported 
the former relationship, while the latter was less certain. Discrimination of translocations from noise 
and other particle activity was addressed by specifically identifying modal translocation peaks in the 
event size distribution, consistently applying a noise threshold across all experiments, and including all 
activity in the event frequency statistics. 

Future experiments will probe the theoretical models and experimental variables more closely, and 
move towards reliable miniaturization of the technology. Improvement of experimental characterization 
procedures will provide the most significant advances towards efficient, accurate pore-sizing. Important 
areas for development include assessment of the internal profile of pores and studies of the charge- 
based solution chemistry, including determination the surface potential on the interior surface of a 
pore. As Henriquez et al. 



23j have suggested, only particle size and concentration can be determined 
from nanopore experiments when geometry and chemistry are poorly defined. Ongoing theoretical 
work relating to blockade sizes should concentrate on incorporating charge effects and pore actuation. 
Further developments relating to geometrical aspects are likely to provide gains in accuracy which are 
insignificant when compared with the experimental spread of values arising from simple models. 
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